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Abstract: In this talk, we consider the interplay of fundamental and matter-induced T-violating 
effects in neutrino oscillations in matter. We present a simple approximative analytical formula for 
the T-violating probability asymmetry for three flavor neutrino oscillations in matter with an arbitrary 
density profile. We also discuss some implications of the obtained results. Since there are no T-violating 
effects in two flavor neutrino case (in the limit of vanshing ^13 or Am|]^, the three flavor neutrino 
oscillations reduces to the two flavor ones), the T-violating probability asymmetry can, in principle, 
provide a way to measure ^13 and Am|]^. 
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1. Introduction 

T violation and CP violation in neutrino oscilla- 
tions have lately been extensively studied in the 
literature [1~3]. However, in most of these stud- 
ies constant matter density was assumed in which 
only fundamental (intrinsic) T violation is feasi- 
ble. To learn more about the effects of CP and 
T violation will be essential and necessary for fu- 
ture experiments such as neutrino factories and 
other long baseline neutrino oscillation experi- 
ments. Future experiments also offer to study 
CPT violation. The measurement of CP, T, and 
CPT violation is very important, not only be- 
cause it will provide us with information about 
neutrino properties, but also because it may have 
interesting implications for physics at high ener- 
gies. 

This talk is based upon the work done by 
E. Akhmedov, P. Huber, M. Lindner, and T. Ohis- 
son [1]. 



* Talk presented at the International Europhysics Con- 
ference on High Energy Physics (EPS HEP 2001), Bu- 
dapest, Hungary, July 12-18, 2001. In collaboration with: 
Evgeny Akhmedov, Patrick Huber, and Manfred Lindner. 
To be published in the Conference Proceedings section of 
the Journal of High Energy Physics (JHEP). 



2. T-violating probability asymme- 
try 

We discuss the interplay of fundamental and matter- 
induced T violation in three flavor neutrino oscil- 
lations in matter.^ T violation cannot be directly 
experimentally tested, since this would mean chang- 
ing the direction of time. However, instead of 
studying neutrino oscillations "backward" in time, 
one can study them forward in time, but with 
the initial and final neutrino fiavors interchanged. 
We define a measure of T violation in neutrino 
oscillations as the following differences, which we 
will call the T-odd probability differences: 



Q/3 



P{y^-^vp)~P{y0^v^), (2.1) 



where P{va vp) is the transition probability 
for Va — > up. Furthermore, we will also denote 
P{i'a Up) by Pa/3- 

For two neutrino flavors there are no T-viola- 
ting effects simply because Pe^ = P^e, which 
means that AP^ = 0. For three neutrino flavors 
the situation is more complicated and we have 
to divide the problem into two separate cases: 
vacuum and matter. In vacuum, we have CPT 



^A T transformation is a time reversal transformation. 
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invariance, which means that we have T viola- 
tion if and only if we have CP violation. In mat- 
ter, the situation is different. Matter is both CP 
and CPT-asymmetric, since it consists of par- 
ticles (nucleons and electrons) and not of their 
antiparticles or, in general, of unequal numbers 
of particles and antiparticles. The matter den- 
sity profiles are, of course, either symmetric or 
asymmetric. Examples of the different types of 
matter density profiles are shown in Figs. |l| and 

I- 

The T-odd probability difference AP^^ (three 
neutrino flavors) has been derived using pertur- 
bation theory to first order in the parameter ^13 
for arbitrary matter density profiles [1] 



AP 



-2S23C233 W {A ~ C*)] 
~ -2S13S23C23 (A - Si2^) 

X 3[e"^^^-r (A, -Ca)] , (2.2) 



where Sab = sin 9 at, Cab = cos 9, 
A 



ab 1 
t 



2_B„ ' 



A = Aa = aJl a*fdt'+Pj^^ (3*fdt', and 

Ca = f J^^^af*dt'. Here a — a{t,to) and /? — 
j3{t,tf)) are to be determined from the solutions 
of the two flavor neutrino problem in the (1, 2)- 
sector (see Ref. [1] for details) and 



/ = f{tM) 
^exp|-z^ (^A-i[5 + y(t')]) dt' 



where V{t) ~ '^^F'^P^^) ^^"^ matter poten- 
tial with Gf ^ 1.16639 • 10~23 eV^^ being the 
Fermi weak coupling constant, tun — 939 MeV 
the nucleon mass, and p{t) the matter density. 
Note that formula (2.2) is valid only when 6*13 
and (5/A are small parameters. In addition, it 
holds in general that AP^^ = APj^ = AP^^ [4], 
i.e., the T-odd probability differences are cyclic 
in the indices and there is in fact only one inde- 
pendent T-odd probability difference. Further- 
more, we have explicitly calculated the T-odd 
probability difference AP^^^ for 

1. matter consisting of two layers (lengths Li 
and L2) of constant density (matter-induced 
potentials Vi and V2) and 

2. an arbitrary matter density profile in the 
adiabatic approximation. 



In the first case in the low energy regime {5 > 
^1,2), we obtain 

sin(26li - 26I2) 



APg^ ~ cos(5cp • 8S12C12S13S23C23 



sin 26'i2 



Jctf 

X {sincjiLi sincj2p2 \Y ~ cos (AiLi + A2P2)]} 

+ sin (5cP • 4si3S23C23 

X Xi [y - cos (AiLi + A2L2)] , (2.3) 

where 9i and 02 arc the matter mixing angles (in 
the (l,2)-sector) in the first and second layers. 



respectively, Aa = A — i (5 + Va) (a 



1,2), 



Y = COSLJiPi COSCiJ2i2 

— sinwiLi sina;2L2 cos(20i — 26*2), 
Xi = sinwiii cosLi;2-^2 sin20i 
-|- sincj2i2 coswiLi sin 2^2 



with 

= ^\/(cOS 2012^ - Vaf + sin^ 2012^2 

(a = 1, 2) and J^s is an effective Jarlskog invari- 
ant similar to the usual Jarlskog invariant [5] 

J = Sl2Cl2Sl3C?3S23C23 sin(5cp. 



In formula (2.2), the cos^cp and sinJcp terms 
describe matter-induced (extrinsic) and funda- 
mental (intrinsic) T violation, respectively. 

In the second case in the regime in which 
the oscillations governed by large A are fast and 
therefore can be averaged over {V{t) < 5 A), 
we obtain 

A r,T c . sin(26l - 26*0) 

AP ~ cos dcp • 4 S12C12S13S23C23 

sm 26'i2 



Jett 

X cos^ $ 

+ sinfcp • 2S13S23C23 sin(6' -f 9^) cos{9 - 9a) 
X sin 2$, (2.4) 

where Oq and 9 are the matter mixing angles (in 
the (l,2)-sector) at the initial and final points 
of the neutrino evolution, to and t, respectively, 
$ = w(t') dt' with 



Ljit) = ^^{cos29uS-Vit)f 



sin^ 26*12(52. 



Note that the cos dcp term can again be written 
in terms of the effective Jarlskog invariant. 
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3. Precision of the analytical approx- 
imation 

Next we discuss the precision of the approxi- 
mative analytical formula for the first case, i.e., 
matter consisting of two layers of constant den- 
sity. In order to do so, we have to distinguish 
between two cases: L/E < 10'' km/GeV and 
L/E > 10^ km/GeV. 

The first case is shown in Fig. ^, whereas the 
second case is shown in Fig. |[ In the first case, 
the oscillating structure of the T-odd probability 
difference AP^ can be resolved. The amplitude 
of APg^ is reproduced very well; however, there 
is an error in the phase that is accumulating with 
the baseline length L as well as with growing 6*13 
and Am^i. In the second case, the oscillations 
governed by the Am|]^ are very fast. The left 
plot uses the same parameter values as in Fig. 
3a of P.M. Fishbane and P. Kaus [2], whereas the 
right plot uses larger values of ^13 and Am^j. 

What about T violation in future terrestrial 
neutrino oscillation experiments including neu- 
trino factory experiments? Will it give rise to 
sizeable effects that are measurable? We studied 
several different experimental setups for a neu- 
trino factory with a beam energy of 50 GeV. Fur- 
thermore, we assumed 2 • 10^^ muon decays and a 
detector mass of 40 kton. (See Ref. [1] for further 
details.) We investigated two different two-layer 
matter density profiles with densities 1 g/cm^, 
3 g/cm'^ and 3 g/cm'^, 3.3 g/cm'^, respectively. 
The first one could correspond to a sea-earth sce- 
nario and the second one to a very long base- 
line experiment in which there should be density 
perturbations. We simulated these scenarios and 
performed fits to the obtained event rates. In ad- 
dition, we compared these simulations for sym- 
metrized versions of the corresponding matter 
density profiles. The symmetrized matter den- 
sity profiles are modeled by replacing the transi- 
tion probabilities with the symmetrized ones: 

- ^(Pdh +Prcv), (3.1) 

where Pdir and Prcv are the transition probabili- 
ties originating from the neutrino propagation in 
the (direct) matter density profile and the "re- 



verse" matter density profile, respectively.^ Thus, 
the simulations are only sensitive to errors in- 
duced by the asymmetry of the matter density 
profile. The difference of the minimal values of 
the functions for the asymmetric and sym- 
metrized matter density profiles is a measure of 
matter-induced T violation. 

Our simulations show that the T-violating 
effects can be "quite sizeable" for the sea-earth 
matter density profile; however, only for L > 
1000 km, which cannot be realized on Earth. The 
qualitative statements of the simulations do not 
change very much if one changes the value of 6c p 
in the fits. For the matter density profile with 
10% density perturbations, the matter-induced 
T-violating effects are small for any baseline. 

4. Summary and conclusions 

In summary, approximative analytical formulas 
for the T-odd probability differences APj^ for 
an arbitrary matter density profile have been de- 
rived using perturbation theory. 

Our main conclusions are the following: 

• T-violating effects can be considered as a 
measure of genuine three-flavorness. 

• For terrestial experiments matter-induced 
T-violating effects can safely be ignored. 

• Asymmetric matter effects cannot hinder 
the determination of the fundamental CP 
and T-violating phase Sep in long baseline 
experiments. 
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Figure 1: Examples of symmetric matter density profiles. Note that also vacuum can be thought of as a 
symmetric matter density profile. 



I <~L 

Figure 2: Example of an asymmetric matter density profile (with two layers of constant density). 
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Figure 3: The T-odd probability difference AP^ for two layers of widths Li = L2 — L/2 and densities 
1 g/cm'^ and 3 g/cm? as a function of the baseline length L. Solid curves show analytical results, whereas 
dashed curves show numerical results. Left plot: On — 0.1 and Am2i = 5 ■ 10^^ eV^. Right plot: On — 0.16 
and Amii = 2 • 10"* eV^. Remaining parameters: Amii = 3.5 ■ 10"^ eV^, 6u = 0.56, 62Z = 7r/4, 5cp = 0, 
and E = 1 GeV. The figure has been taken from Ref. [1]. 
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Figure 4: The same as in Fig. ^ except that the densities are and 6.4 g/cm'^ and E — 0.5 GeV. The grey 
curves are analytical results, whereas the black solid and dashed curves show the results averaged over fast 
oscillations of the analytical and numerical calculations, respectively. The figure has been taken from Ref. [1]. 
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